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ABSTRACT 


A global  model  is  presented  which  predicts  the  mass  burning  flux  for  pool  fires  consuming 
liquid  fuels  in  a quiescent  environment.  The  model  assumes  constant  bulk  properties  such  as  flame 
temperature,  soot  volume  fraction,  and  species  concentration.  The  computational  procedure 
requires  knowledge  of  the  fuel  smoke  point  height  and  fuel  properties  such  as  the  heat  of 
vaporization,  heat  capacity,  and  boiling  point.  A cylindrical  flame  shape  is  assumed  with  the  flame 
height  given  by  Heskestad's  correlation.  A mean  beam  length  approach  for  radiative  heat  transfer  is 
utilized  and  emission  from  both  gas  species  and  soot  particles  is  considered.  The  convective  heat 
transfer  coefficient  is  estimated  using  a Raleigh  number  correlation.  Experiments  in  small  diameter 
pool  fire  are  used  to  quantify  the  conductive  heat  transfer.  The  predicted  mass  flux  for  a number  of 
fuels  is  within  a factor  of  two  of  measured  burning  rates  for  pool  diameters  greater  than  0.2  m. 
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INTRODUCTION 


The  mass  flux  of  burning  fuel  in  a pool  fire  depends  on  the  heat  feedback  from  the  flame  to 
the  fuel  surface,  which  occurs  in  the  form  of  conductive,  convective,  and  radiative  heat  transfer. 
The  relative  contribution  of  each  of  these  components  depend  on  a large  number  of  parameters 
including  pool  diameter,  burner  material,  lip  height,  flame  shape,  and  the  spatial  distribution  of 
temperature,  species  concentration  and  soot  volume  fraction.  The  goal  of  the  present  study  is  to 
develop  a predictive  algorithm  for  burning  flux  which  is  consistent  with  the  essential  physics  of  the 
problem  and  which  is  also  generic,  such  that  experimental  determination  of  parameters  in  the  model 
is  not  required  for  each  new  fuel  or  pool  diameter  considered. 

A number  of  burning  rate  models  have  been  developed  using  a global  approach. 
[1,2,3 ,4,5, 6]  The  heat  feedback,  Q (W),  to  the  pool  surface  is  the  sum  of  conduction  ( Qcond ), 
convection  ( Qcom  ),  and  radiation  ( ): 


Q Qcond  Qconv  Qrad 

(1) 

The  model  considers  the  fuel  to  be  isothermal  horizontally  and  to  be  steadily  burning.  Hottel  [1] 
represented  heat  transfer  to  the  pool  surface  in  terms  of  characteristic  or  global  flame  properties: 

Qcond  = k ft  D(Tf  -Ts) 

(2) 

Qconv  =hAs(Tf  - Ts ) 

(3) 

Qrad  =gF As(Tj  - e~rD) 

(4) 

where  k (W/m-K)  is  a conduction  coefficient,  D (m)  is  the  pool  diameter. 

As  (m2)  is  the  pool 

surface  area,  Tf  (K)  is  the  flame  temperature,  Ts  (K)  is  the  pool  surface  temperature,  h (W/m2  K)  is 

-8  ? 

a convective  heat  transfer  coefficient,  a is  the  Stefan-Boltzmann  constant  (5.67 x 1 0 W/m  K),  F is 
a dimensionless  flame-pool  surface  radiative  view  factor,  and  p (m'1)  is  a radiative  extinction 
coefficient. 

Hottel  noted  that  when  D is  small,  conduction  dominates  because  convection  and  radiation 
are  proportional  to  D 2 . When  D is  large  conduction  becomes  unimportant  and  radiation  eventually 
dominates  convection,  because  pD  becomes  large  and  radiation  is  proportional  to  T* . 

Global  burning  rate  models  have  been  discussed  by  a number  of  other  researchers. 
Hayasaka  and  Koseki  [2]  used  a mean  beam  length  model  with  global  flame  properties  and 
compared  their  predictions  with  measured  burning  rates  for  pools  up  to  10  m in  diameter.  Their 
model,  however,  has  not  been  tested  on  fuels  other  than  kerosene.  Modak  and  Croce  [3]  calculated 
the  radiative  heat  flux  from  the  flame  to  a surface  element  on  a burning  pool  (0.18  m)  of  PMMA  in 
terms  of  an  empirical  time-averaged  flame  shape,  an  effective  radiation  temperature,  and  a mean 
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gray-body  absorption-emission  coefficient.  These  parameters  were  all  independently  measured 
whereas  convection  was  assumed  to  be  negligible.  Orloff  and  de  Ris  [4,5]  represented  the  fire  as  a 
time  averaged  volume  of  constant  property  gases.  Their  measured  radiative  heat  loss  fraction  (Xr) 
was  used  to  estimate  a global  absorption-emission  coefficient.  Measured  flame  shapes  were  used  to 
calculate  a mean  beam  length.  Convective  heat  feedback  was  estimated  from  experiments  on  a 
water  cooled  gas  burner.  Model  predictions  were  compared  to  measurements  for  several 
experimental  conditions  for  pools  burning  liquid  methanol  and  solid  thermoplastics.  An  attempt  to 
validate  their  model  for  liquid  fuels  other  than  methanol  was  not  successful.  This  may  be  because 
flame  shapes  can  vary  dramatically  for  different  fuels  (as  a function  of  stoichiometry)  and  thus  time 
averaged  shapes  must  be  determined  for  each  fuel  of  interest. 

MODEL  FORMULATION 

We  consider  the  mass  burning  rate  of  a liquid  fuel  in  a metal  burner  of  internal  diameter  D. 
The  case  of  small  lip  height  (/)  and  moderate  pool  depth  (d)  is  examined  (see  Fig.  1).  Heat  flux 
(W/m2)  feedback  to  the  pool  surface  is  based  on  Eq.  1.  The  burning  rate  is  taken  as  the  sum  of 
conduction  ( Qcond ),  convection  (Qconv),  radiation  (Q^),  and  re-radiation  {Qnrad)  from  the  pool 
surface: 


M[HV  + Cpf  (Ts  -Ta)]~  Qcond  + Qcom  + Qrad  - Qrerad  ( 5 ) 

The  heat  feedback  represents  the  energy  required  to  vaporize  the  liquid  fuel  where  M (g/s)  is  the 
mass  burning  rate,  Hv  (J/g)  is  the  heat  of  vaporization  at  Ts  (K),  Cpf  (J/g  K)  is  the  liquid  heat 
capacity,  and  Ta  (K)  is  the  bulk  fluid  temperature  at  ambient  condition.  Equation  5 assumes  that  all 
heat  to  the  pool  surface  goes  to  vaporize  the  fuel,  i.e.,  heat  losses  for  these  deep  pools  (<i>  10  cm) 
and  increases  in  the  sensible  heat  of  the  liquid  fuel  are  negligible.  In  our  model,  Qcond  represents 
heat  transfer  from  the  fire  through  the  burner  walls  into  the  liquid  pool. 

The  flame  is  modeled  as  a homogeneous  mixture  of  entrained  air,  combustion  products 
(C02,  H20,  and  soot),  unbumt  fuel,  and  inert  N2  at  a global  flame  temperature.  The  global 
temperature,  Tf  , is  calculated  from  an  enthalpy  balance  about  a control  volume  encompassing  the 

gas-phase  cylindrical  flame.  The  amount  of  entrained  air  entering  the  flame  is  determined  from 
Reference  [7],  modified  by  a constant  factor  which  is  determined  from  a best  fit  of  the  burning  rate 
data.  This  form  of  the  global  temperature  was  chosen  only  after  trying  a number  of  approaches  to 
find  a suitable  representative  global  temperature,  including  a constant  temperature  (1000  K,  1100 
K,  1200  K,  etc  ) and  a temperature  based  on  the  average  of  the  calculated  maximum  flame 
temperature  and  the  surface  temperature  (7^),  which  was  assumed  to  be  at  the  boiling  point.  These 
forms  led  to  poor  predictions  of  the  burning  rate. 

Tj  is  calculated  in  a manner  similar  to  Smith  and  Van  Ness  [8].  The  effects  of  radiative 
heat  losses,  combustion  efficiency,  and  dilution  by  air  entrainment  are  also  considered.  The  energy 
available  to  heat  the  excess  entrained  air,  combustion  products,  and  unbumt  fuel  to  7)  is  equal  to 

the  heat  of  combustion,  Hc  (J/mol),  less  the  energy  due  to  losses  associated  with  radiative  loss 
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(ZR  ),  combustion  efficiency  (Za ),  and  the  heat  of  vaporization  of  the  fuel,  Hv  (J/mol),  at  ambient 
temperature,  Ta  (K): 


X>Cp)(T)dT  (6) 

/ ° 

where  Cp  i (7)  (J/mol  K)  is  the  vapor  heat  capacity  [9]  of  the  F gaseous  species,  which  includes 
N2,  C02,  H20,  unbumt  fuel,  and  excess  02,  and  X\  is  the  mole  fraction  of  the  gaseous  species. 
The  radiative  heat  loss  fraction  is  defined  as: 


y - ( 7 \ 

R MHC  / MW 

where  Qr  (W)  is  the  energy  emitted  by  the  flame  to  the  surroundings.  The  combustion  efficiency  is 
defined  as: 


// 


MH/MW 


(8) 


where  Ha  (W)  is  the  actual  heat  release  rate  and  M ■ Hc  / MW  (W)  represents  the  ideal  heat 
release,  where  MW  is  the  fuel  molecular  weight  (g/mol).  The  concentrations  ( Xi ) of  C02,  H20,  N2 
and  02  are  determined  from  stoichiometry  and  the  amount  of  entrained  air.  The  concentration  of 
unbumt  fuel  is  determined  from  the  combustion  efficiency.  The  chemical  reaction  for  1 mole  of 
fuel  (F)  containing  V carbon  atoms,  "h"  hydrogen  atoms,  and  "o"  oxygen  atoms  mixing  with 
entrained  air  in  the  flame  is  represented  as: 


F + a[0.2\O2  + 0.79  /V,]  — » 


c CO 2 +-H20 
2 2 2 


+ (1  - Xa)F  + 0.79  aN2+  (0.21  a - kZa)02 


(9) 


where  A:  is  a constant  = [c  + (h/4)  - (o/2)]  and  a is  the  moles  of  air  (=  Ma/[28.96  g/mol])  entrained 
into  the  fire.  The  mass  of  air  per  unit  time  ( Ma  ) entrained  into  a fire  of  height  H,  is  given  by  [7] 
(in  units  of  kg/s): 


Ma  = 0.08 6(H)05  (rZa  + 1 )D2  Frc' 

for  H1  < 1 

(10) 

Ma  = 0.093  (if)1 5 ( rZa  +1  )D2  Frc' 

for  1 < R < 4 

(11) 

Ma  =0.018(/T)25  (rZa  + \)D2  Frc' 

for  77  > 4 

(12) 
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where  H'  (=  HID)  is  a non-dimensional  flame  height  and  c'  is  an  empirically  determined 
dimensionless  parameter  (=  0.11),  universal  for  all  fuels,  and  r is  the  stoichiometric  mass  ratio  of 
fiile  to  air.  The  Froude  number  is  defined  as: 


Fr  = 


QpR' 

\Hc(Xo-Xr) 

V R' 


where  R'  = (r  xa  + 1)  Cpa  Ta  and  Q'D  = 


MHC 

p.  c r.  D 2 


(13) 


Intermediate  gas  species  (CO,  H2,  etc.)  and  soot  are  neglected  in  the  temperature 
calculation,  except  in  a global  sense  through  xa  • Table  1 summarizes  some  of  the  fuel  parameters 
used  in  the  temperature  calculation  including  the  fuel  type,  the  radiative  heat  loss  fraction  (%R) 
[6,10,11],  the  combustion  efficiency  (xa)  [10]  and  the  characteristic  soot  volume  fraction  (/„) 
which  is  discussed  below.  The  calculated  flame  temperature  ( Tf ) is  a weak  function  of  diameter  as 
the  fuel/air  entrainment  varied.  Tf  for  a 1 m diameter  pool  is  given  in  Table  1 . 


Table  1 Summary  of  fuel  parameters  used  in  the  model  including  fuel  type,  radiative  heat 
loss  fraction  (zx),  combustion  efficiency  ( j0),  soot  volume  fraction  (/v),  and  the 
calculated  global  flame  temperature  ( Tf  ) for  aim  diameter  pool  fire. 


Fuel 

** 

fv  x 10^ 

Tf  (K) 

Methanol 

0.20b 

1.0a 

0.0 

1114 

Ethanol 

0.23a 

0.92a 

0.07 

1034 

Acetone 

0.26a 

0.94a 

0.11 

1042 

MMA 

0.34b 

0.88a 

0.22 

933 

Hexane 

0.30a 

0.93a 

0.29 

1007 

Heptane 

0.30b 

0.92a 

0.50 

1010 

Toluene 

0.32b 

0.76d 

6.0 

843 

Benzene 

0.35c 

0.76e 

6.8 

802 

a Reference  [10] 
b Reference  [11] 
c Reference  [6] 

d measured  in  cone  calorimeter  [12] 
e assumed  equal  to  xa  f°r  toluene 
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Radiation 


Heat  transfer  to  the  surface  of  the  pool  by  radiation  is  modeled  as  isothermal  emission  from 
gaseous  C02  and  H20  and  gray-body  emission  from  soot.  A correction  for  spectral  overlap  by 
these  emitters  was  included.  Total  radiation  to  the  pool  surface  ( Qrad ) was  given  by  Siegel  and 
Howell  [13]: 


Qrad  = cjT;  As  (1  - p)[  1 - (1  - esoot)(\  - 0(1 " *c)]  ( 14  ) 

where  p (=  0.08)  is  the  surface  reflectivity  for  liquid  fuels  [14],  £soot  is  the  spectrally  averaged 
emissivity  of  soot,  sw  is  the  total  emissivity  of  H20,  and  sc  is  the  total  emissivity  of  C02.  Gas 
emission  and  absorption  by  fuel  vapor,  CO  and  other  gas  species  is  not  considered  here.  Gas 
emissivity  for  C02  and  H20  are  calculated  using  the  algorithm  given  by  Leckner  [15].  Reflection 
of  the  incoming  radiation  is  dependent  on  the  angle  of  incidence.  A value  of  0.08  is  a reasonable 
value  for  pool  fires,  regardless  of  fuel  type  [14]. 

Yuen  and  Tien  [16]  have  shown  that  soot  emission  can  be  approximated  by  the  following 
expression  for  a gray  emitter: 


£soo,  =l-e 


-K  L 


(15) 


In  the  above  equation,  L is  the  path  length  and  k is  the  effective  soot  emission  parameter: 


k = 


3 .6(3,2} 

Q 


(16) 


where  C2  is  Plank's  second  constant  (0.014388  m K)  and  the  dimensionless  effective  soot 
concentration  parameter,  C„  is  defined  as: 

36tc  fv  n2  s 

Q =r ^ (17) 

[«2  - (ns)2  +2]  + 4 n4  s2 

The  infrared  averaged  optical  constants  (n  = 3.49  and  s = 2.17)  for  soot  particles  were  taken  from 
Tien  and  Lee  [17].  For  a particular  flame  geometry,  the  assumption  of  a homogeneous  gray 
mixture  allows  use  of  the  mean  beam  length  ( Lm ) to  represent  the  characteristic  path  length,  L,  (see 
Eq.  15).  For  simplicity,  the  flame  is  assumed  to  be  a right  cylinder,  and  the  mean  beam  length  is 
given  by  [5]: 


L = 


3.6  V 


(18) 
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where  V is  the  flame  volume  and  Af  is  the  flame  surface  area.  For  a right  cylinder,  the  flame 
volume  and  surface  area  are  given  by: 


V = tcR2  H (19) 

Aj  = 2nR(R  + H)  (20) 

where  R is  the  pool  radius  and  H is  the  flame  height  determined  using  the  correlation  given  by 
Heskestad  [18]. 

A number  of  studies  have  measured  soot  volume  fraction  distributions  in  pool  flames 
burning  liquid  fuels  [11,19]  whereas  only  a few  studies  have  systematically  studied  the  effects  of 
scale  or  fuel  type.  Bard  and  Pagni  [20]  measured  the  soot  volume  fraction  0.03  m above  the  surface 
of  0.20  m diameter  pool  fires  for  a number  of  fuels.  Whereas  /v  changes  as  a function  of  location 
in  the  flame,  Bard  and  Pagni's  measurements  were  taken  as  representative  of  the  soot  volume 
fraction.  For  fuels  that  Pagni  did  not  measure,  a correlation  with  smoke  point  height  (ls)  is  used. 
The  smoke  point  height  characterizes  the  sooting  propensity  of  fuels  and  has  been  measured  for  an 
extensive  number  of  fuels  [10,21,22].  Figure  2 shows  Pagni's  /„  measurements  as  a function  of 
fuel  smoke  point  height.  Although  it  is  possible  that  combustion  characteristics  may  vary  as  a 
function  of  pool  diameter,  not  enough  reliable  measurements  exist  to  determine  that  dependence. 
The  existing  data  suggest  that  for  pool  diameters  less  than  2 m,  xa  and  Xr  are  independent  of  scale 
[23].  Because  of  the  relationship  between  /v  and  xa  > 35  shown  in  Fig.  2,  the  model  also  assumes 
that  /„  is  a function  of  fuel  type  only.  Variation  of  combustion  characteristics  as  a function  of  pool 
diameter  can  be  easily  incorporated  into  the  present  model. 

Because  of  the  relatively  low  boiling  point  of  these  fuels,  re-radiation  from  the  pool  surface 
to  the  ambient  [4,5]: 


Qrerad=°MT?  -T?)  (21) 

is  very  small  compared  to  convection  or  radiation. 

Convection 

Convective  heat  transfer  to  the  fuel  surface  is  approximated  using  the  stagnant  film  model 

[4,5]: 


Qconv 


= A . 


h' 


C 


pa 


Hc  (Xa  ~ Xr)  r 

Xa 


-Cpa{Ts-Ta) 


(22) 


where  y = M(Cpa/hy  As  is  a blowing  factor  which  effectively  reduces  heat  transfer  as  mass 
transfer  increases.  The  symbol  r represents  the  stoichiometric  air/fuel  mass  ratio  and  Cpa  is  the 
effective  heat  capacity  of  air  taken  at  a representative  temperature  equal  to  ( 7}  + Ts)/ 2.  A simple 
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heat  transfer  coefficient  (h')  is  used,  based  on  a Rayleigh  number  criterion  for  natural  convective 
heat  transfer  on  a hot  surface  facing  upwards  [24]: 

h'  = 054Ra025  \05<Ra<2-\07 

h'  = 0.14ita°'33  2-107  < Ra  < 3-1010  ( 23  * 

h'  takes  on  values  typically  « 9 W/m2  K,  which  is  consistent  with  the  value  used  by  Orloff  and  de 
Ris  (h'&  7 W/m2  K)  for  pools  with  finite  height  lips  [4,5].  If  convective  transfer  to  the  pool  surface 
is  modeled  as  flow  past  a cavity,  then  the  effect  of  lip  height  on  h'  is  small  [25],  and  use  of  the  same 
value  of  h'\o  model  convective  transfer  in  a fire  without  a finite  lip  height  is  not  unreasonable. 

Conduction 

Conduction  from  the  flame  to  the  fuel  surface  is  modeled  as  one-dimensional  heat  transfer 
vertically  down  through  the  burner  walls.  All  heat  is  assumed  to  be  transferred  into  the  liquid  pool: 

dT 

Qcond  - kb  1 ( 24  ) 

d z 


where  kb  (W/m  K)  is  the  thermal  conductivity  of  the  burner  and  is  material  dependent,  d77dz 
(K/m)  is  the  thermal  gradient  in  the  wall  and  Aw  (m2)  is  the  wall  cross  sectional  area: 


. *(■ D2o~D 2) 

- 4 

where  Da  (m)  is  the  outer  diameter  and  D (m)  the  inner  diameter  of  the  pan. 


(25) 


To  determine  the  thermal  gradient  dT/dz  in  Eq.  24,  temperature  measurements  were  made 
inside  the  walls  of  brass  burners.  The  experimental  set-up  has  been  described  previously  [11]. 
Thermocouples  were  spaced  0.003  m apart  (0.001  m,  0.004  m,  and  0.007  m below  the  burner  rim) 
and  in  the  center  of  the  brass  wall  of  a 0.035  m diameter  pool.  During  each  bum  a constant  lip 
height  (0.007  m)  was  maintained.  A temperature  gradient  of  160  K/m  was  measured  a few  seconds 
after  ignition  which  remained  constant  during  the  entire  experiment  (20  min).  Decreasing  the  lip 
height  to  0.005  m did  not  change  the  temperature  gradient.  Nor  did  the  temperature  gradient 
change  for  a 0.050  m diameter  pool  with  0.005  m lip  height  and  0.0032  m brass  walls. 
Thermocouples  mounted  0.001  m and  0.002  m into  the  walls  0.001  m below  the  burner  rim 
measured  the  same  temperatures  (within  0.1  K)  as  a function  of  time  from  flame  ignition,  indicating 
that  the  assumption  of  one-dimensional  vertical  heat  conduction  is  appropriate.  Table  2 lists  the 
conditions  and  results  for  ethanol  burning  rate  experiments.  The  total  conductive  heat  transfer 
down  through  the  burner  walls  increased  for  the  larger  pool,  but  the  conductive  flux  normalized  by 
the  pool  surface  area  increased  only  slightly. 
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Using  Eq.  24,  the  total  energy  conducted  ( Qcond  ) to  the  fuel  increases  with  pool  diameter 
(D)  as  noted  in  Table  2,  whereas  the  energy  conducted  per  unit  wall  area  ( Q"  ) is  approximately 
constant. 


Table  2 Pool  diameter  (D),  burner  wall  thickness,  measured  thermal  gradient  inside  of  the 
burner  walls,  total  conductive  heat  transfer  through  the  burner  walls  ( Qcond  ),  and 

average  conductive  heat  transfer  per  unit  area  of  pool  surface  ( Q " ) for  pools  with 
brass  walls  burning  ethanol. 


D{  m) 

Wall  (m) 

dT/dz  (K/m) 

fi-  (W) 

Q"  (W/m2) 

0.035 

0.0020 

160  ±10* 

3.9 

4100  ±400 

0.050 

0.0032 

160  + 10 

8.8 

4500  ±500 

*the  uncertainty  is  expressed  as  1 a. 

RESULTS  AND  DISCUSSION 

The  required  input  for  the  model  are  physical  parameters  such  as  the  fuel  heat  of 
vaporization,  liquid  and  vapor  fuel  heat  capacity,  the  fuel  boiling  temperature,  the  burner  material 
composition  and  wall  thickness,  and  the  combustion  characteristics  of  the  fuel  such  as  the 
stoichiometry,  combustion  efficiency,  the  radiative  heat  loss  fraction  and  the  soot  volume  fraction. 
The  model  was  tested  on  a variety  of  liquid  fuels  over  a range  of  pool  diameters  from  0.1  m to  2.5 
m.  The  fuels  tested  include  methanol,  ethanol,  acetone,  methyl  methacrylate  (MMA),  hexane, 
heptane,  benzene  and  toluene.  These  fuels  yield  flames  with  a wide  range  of  flame  luminosities, 
heat  release  rates,  and  soot  volume  fractions. 

The  predicted  mass  burning  flux  was  calculated  by  solving  Eq.  5 iteratively.  The  range  of 
pool  diameters  considered  was  from  0.1  to  2.5  m.  The  upper  bound  was  limited  by  measurements 
available  in  the  literature.  Conduction  was  modeled  by  assuming  that  the  burner  material  was  brass 
with  walls  0.002  m thick.  These  conditions  were  chosen  as  representative  in  order  to  quantitatively 
model  conduction,  and  realizing  that  the  conditions  vary  widely  from  experiment  to  experiment. 
The  calculations  showed  that  conduction  is  small  compared  to  convection  and  radiation. 

Table  3 lists  values  of  the  ratio  of  the  mean  beam  length  to  pool  diameter,  the  ratio  of  flame 
height  to  pool  diameter,  and  the  emissivities  of  soot,  C02  and  H20  as  a function  of  fuel  type  and 
pool  diameter.  For  each  fuel,  the  ratio  of  mean  beam  length  to  pool  diameter  is  nearly  independent 
of  diameter  indicating  that  the  absolute  value  of  the  mean  beam  length  is  proportional  to  pool 
diameter.  The  ratio  of  flame  height  to  diameter  decreases  weakly  as  a function  of  pool  diameter. 
The  flame  heights  and  the  mean  beam  lengths  for  methanol  are  much  smaller  than  heptane  or 
toluene.  The  values  of  the  emissivities  are  functions  of  stoichiometry,  temperature  and  mean  beam 
length  and  thus  vary  with  fuel  and  pool  diameter.  As  expected,  the  soot  emissivity  for  toluene  is 
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the  highest  whereas  the  gaseous  emissivities  for  the  three  fuels  are  comparable.  The  magnitude  of 
the  values  are  in  line  with  those  determined  previously  [26]. 

Figures  3-10  compare  measured  mass  burning  fluxes  to  those  predicted  by  the  model  as  a 
function  of  pool  diameter  for  the  fuels  listed  in  Table  1 . The  experimental  measurements  were 
compiled  from  a large  number  of  studies  [1 1, 14, 19„27,28, 29,30,3 1,32].  The  scatter  in  the 
experimental  data,  as  in  Fig.  8 for  example,  may  be  attributed  to  differences  in  the  experimental 
method,  burner  material,  environmental  conditions  (quiescent  or  windy),  lip  height,  use  of  water 
cooling,  etc.  Whereas  a few  experiments  have  been  performed  under  quasi  "steady  state" 
conditions  with  a constant  lip  height  accomplished  by  replenishing  the  fuel  (filled  symbols  in  Figs. 
3-10),  most  experiments  have  been  conducted  under  free  bum  conditions  where  the  lip  height 
increases  with  time  (open  symbols  in  Figs.  3-10).  As  the  lip  height  increases,  large  changes  in  the 
radiative  heat  transfer  can  occur  as  the  view  angle  from  flame  to  fuel  surface  decreases.  The 
convective  heat  transfer  may  also  decrease.  Because  the  number  of  constant  lip  height  burning  flux 
measurements  in  the  literature  are  limited,  additional  measurements  were  made  using  the  apparatus 
and  methods  described  previously  [1 1].  In  Figs.  3-10,  both  constant  and  varying  lip  height  data  are 
plotted.  Steady  Up  height  burning  rate  data  are  indicated  by  the  filled  symbols  and  bum-down 
(varying  lip  height)  data  are  denoted  by  the  hollow  symbols.  Only  a few  measurements  are 
available  to  compare  burning  fluxes  under  conditions  of  constant  versus  varying  lip  height  (for  the 
same  diameter  pool).  Fig.  3 (D  = 0.3  m and  0.6  m),  Fig.  4 (D  « 0.3  m),  Fig.  8 (D  = 0.6  m)  and  Fig. 
10  (D  = 0.6  m)  show  the  same  burning  flux  regardless  of  lip  height  conditions  whereas  only  Fig.  8 
( D = 0.3  m)  shows  some  difference  in  burning  fluxes.  For  these  reasons,  all  available 
measurements  are  compared  to  the  calculations  presented  here. 

Figures  3-10  show  that  the  burning  rate  model  calculations  agree  with  experiment  to  within 
approximately  a factor  of  2 for  pool  diameters  greater  than  0.2  m.  The  largest  differences  occur  for 
small  pool  diameters  (<  0.2  m)  where  the  model  was  within  a factor  of  three  or  better. 

The  sensitivity  of  the  burning  flux  calculation  to  variation  of  input  parameters  was  assessed 
by  defining  a sensitivity  coefficient  (Sf)  such  that: 


A ml  m 

Apt/ Pi 


(26) 


where  (A m/m)  was  the  percentage  change  in  the  mass  burning  flux  and  (A  p,/p,)  was  the 
percentage  change  in  one  of  the  input  parameters  %a->  Xr>  /,>or  Tf  . The  sensitivity  coefficients 
determined  for  the  0.30  m heptane  and  methanol  pools  are  presented  in  Table  4.  S,  varies  only 
slightly  with  pool  diameter  for  each  of  these  fuels.  A negative  sensitivity  coefficient  indicates  a 
decrease  in  mass  burning  flux  for  an  increase  in  the  perturbed  parameter.  Table  4 shows  that  the 
largest  change  in  burning  rate  occurs  for  a perturbation  of  temperature,  followed  by  ja,  %R,  and 
/v.  If  the  global  temperature  estimate  is  good  to  150  K,  then  the  burning  flux  calculation  is 
accurate  to  50  % for  heptane  and  40  % for  methanol. 

Figure  11  compares  the  calculated  and  measured  [14]  percentage  of  radiative  transfer  as  a 
function  of  pool  diameter  for  methanol,  heptane  and  toluene  fires.  For  all  fuels,  the  relative 
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contribution  of  radiation  increases  with  pool  diameter.  As  expected  from  the  %R  data,  radiative 
heat  transfer  to  the  pool  surface  is  most  important  for  toluene  followed  by  heptane  and  methanol 
respectively,  for  all  pool  diameters. 


Table  3 The  calculated  values  of  the  ratio  of  the  mean  beam  length  ( Lm)  to  pool  diameter 
(Z>),  the  ratio  of  flame  height  ( H)  to  pool  diameter  ( D ),  and  the  emissivities  of  soot 
( £soc  )>  HzO  ( £w  ) and  C02  ( sc ). 


Fuel 

D(  m) 

LJD 

HID 

p 

° soot 

£c 

Methanol 

0.05 

0.72 

2.0 

0 

0.028 

0.037 

0.1 

0.70 

1.8 

0 

0.046 

0.049 

0.5 

0.68 

1.5 

0 

0.13 

0.081 

1 

0.66 

1.4 

0 

0.18 

0.097 

Heptane 

0.05 

0.82 

4.8 

0.014 

0.023 

0.042 

0.1 

0.81 

4.3 

0.028 

0.039 

0.055 

0.5 

0.80 

4.2 

0.13 

0.11 

0.088 

1 

0.80 

4.1 

0.25 

0.16 

0.10 

Toluene 

0.05 

0.80 

4.2 

0.13 

0.020 

0.047 

0.1 

0.80 

4.1 

0.24 

0.034 

0.060 

0.5 

0.80 

4.0 

0.75 

0.095 

0.093 

1 

0.79 

3.7 

0.94 

0.14 

0.11 

Table  4 Sensitivity  of  the  Model  to  Input  Parameters  for  0.3  m diameter  pool  fires. 


S' 

Parameter 

Heptane 

Methanol 

Zr 

-0.90 

-0.5 

Za 

2.1 

1.8 

/v 

0.2 

* 

Tf 

3.5 

2.9 

* S,  not  defiled  because  fv=  0 for  methanol  pool  fires. 
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CONCLUSIONS 


A simple  global  model  for  predicting  the  burning  rates  of  liquid  fuels  has  been  developed. 
The  model  requires  the  input  of  fuel  properties  (Cp/ , Hv,  boiling  temperature,  etc.)  and  fuel 

combustion  characteristics  {x0->Xr^  and  /v ).  Smoke  point  height  measurements  for  an  extensive 
number  of  fuels  are  available  in  the  literature  as  are  correlations  of  xa>  Xr>  and  /„  with  h 
[10,21,23].  For  other  fuels,  a single  simple  experiment  would  be  necessary  to  obtain  ls  [33].  Thus, 
the  model  described  here  requires  no  input  other  than  physical  properties  of  the  fuel  and  the  fuel 
smoke  point  height. 

Although  a comparison  of  the  results  of  the  model  with  measurements  are  encouraging, 
further  work  is  required  to  systematically  measure  the  mass  burning  flux  for  large  (>  0.60  m) 
diameter  pools  under  conditions  of  steady  lip  heights  where  few  or  no  data  currently  exist.  In 
addition,  the  importance  of  radiative  absorption  by  fuel  vapor  is  under  study  by  use  of  absorption 
spectroscopy. 
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Figure  1.  Schematics  showing  various  heat  transfer  mechanisms  associated  with  a liquid  pool 
fire. 
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Figure  2.  Characteristic  value  of  the  soot  volume  fraction  in  a pool  fire  [20]  as  a function  of 
the  smoke  point  height  [10]  for  a number  of  different  fuels. 
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Figure  3.  Comparison  of  experimentally  measured  mass  burning  fluxes  to  those  predicted  by 
the  model  for  methanol  as  a function  of  pool  diameter.  Filled  symbols  represent 
measurements  where  the  lip  height  was  maintained  at  a constant  value.  Open  symbols 
represent  measurements  where  the  lip  height  was  varying. 
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Figure  4.  Comparison  of  experimentally  measured  mass  burning  fluxes  to  those  predicted  by 
the  model  for  ethanol  as  a function  of  pool  diameter.  Filled  symbols  represent  measurements 
where  the  lip  height  was  maintained  at  a constant  value.  Open  symbols  represent 
measurements  where  the  lip  height  was  varying. 
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Figure  5.  Comparison  of  experimentally  measured  mass  burning  fluxes  to  those  predicted  by 
the  model  for  acetone  as  a function  of  pool  diameter.  Filled  symbols  represent  measurements 
where  the  lip  height  was  maintained  at  a constant  value.  Open  symbols  represent 
measurements  where  the  lip  height  was  varying. 
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Figure  6.  Comparison  of  experimentally  measured  mass  burning  fluxes  to  those  predicted  by 
the  model  for  MMA  as  a function  of  pool  diameter.  Filled  symbols  represent  measurements 
where  the  lip  height  was  maintained  at  a constant  value. 
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Figure  7.  Comparison  of  experimentally  measured  burning  rates  to  those  predicted  by  the 
model  for  hexane  as  a function  of  pool  diameter.  Filled  symbols  represent  measurements 
where  the  lip  height  was  maintained  at  a constant  value.  Open  symbols  represent 
measurements  where  the  lip  height  was  varying. 
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Figure  8.  Comparison  of  experimentally  measured  burning  rates  to  those  predicted  by  the 
model  for  heptane  as  a function  of  pool  diameter.  Filled  symbols  represent  measurements 
where  the  lip  height  was  maintained  at  a constant  value.  Open  symbols  represent 
measurements  where  the  lip  height  was  varying. 
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Figure  9.  Comparison  of  experimentally  measured  mass  burning  fluxes  to  those  predicted  by 
the  model  for  benzene  as  a function  of  pool  diameter.  Open  symbols  represent  measurements 
where  the  lip  height  was  varying. 
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Figure  10.  Comparison  of  experimentally  measured  mass  burning  fluxes  to  those  predicted 
by  the  model  for  toluene  as  a function  of  pool  diameter.  Filled  symbols  represent 
measurements  where  the  lip  height  was  maintained  at  a constant  value. 
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Figure  11  Comparison  of  the  measured  [14]  and  calculated  percentage  of  heat  feedback 
which  is  radiative  as  a function  of  pool  diameter  for  methanol,  heptane,  and  toluene  pool 
fires. 
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A global  model  is  presented  which  predicts  the  mass  burning  flux  for  pool  fires  consuming  liquid  fuels  in  a quiescent  environment.  The 
model  assumes  constant  bulk  properties  such  as  flame  temperature,  soot  volume  fraction,  and  species  concentration.  The  computational 
procedure  requires  knowledge  of  the  fuel  smoke  point  height  and  fuel  properties  such  as  the  heat  of  vaporization,  heat  capacity,  and  boiling 
point.  A cylindrical  flame  shape  is  assumed  with  the  flame  height  given  by  Heskestad's  correlation.  A mean  beam  length  approach  for 
radiative  heat  transfer  is  utilized  and  emission  from  both  gas  species  and  soot  particles  is  considered.  The  convective  heat  transfer  coefficient 
is  estimated  using  a Raleigh  number  correlation.  Experiments  in  small  diameter  pool  fire  are  used  to  quantify  the  conductive  heat  transfer. 
The  predicted  mass  flux  for  a number  of  fuels  is  within  a factor  of  two  of  measured  burning  rates  for  pool  diameters  greater  than  0.2  m. 
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